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ABSTRACT 

We present 3D hydrodynamical simulations of ram pressure stripping of a disc galaxy 
orbiting in a galaxy cluster. In this paper, we focus on the properties of the galaxies' 
tails of stripped gas. The galactic wakes show a flaring width, where the flaring angle 
depends on the gas disc's cross-section with respect to the galaxy's direction of motion. 
The velocity in the wakes shows a significant turbulent component of a few 100 kms^^. 
The stripped gas is deposited in the cluster rather locally, i.e within ^ 150 kpc from 
where it was stripped. We demonstrate that the most important quantity governing 
the tail density, length and gas mass distribution along the orbit is the galaxy's mass 
loss per orbital length. This in turn depends on the ram pressure as well as the galaxy's 
orbital velocity. 

For a sensitivity limit of 10^^ cm"^ in projected gas density, we find typical 
tail lengths of 40 kpc. Such long tails are seen even at large distances (0.5 to 1 Mpc) 
from the cluster centre. At this sensitivity limit, the tails show little flaring, but a 
width similar to the gas disc's size. 

Morp hologically, we find g ood agreement with the HI tails observed in the Virgo 
cluster bv IChung et al.l (j2007 ). However, the observed tails show a much smaller ve- 
locity width than predicted from the simulation. The few known X-ray and Ha tails 
are generally much narrower and much straighter than the tails in our simulations. 
Thus, additional physics like a viscous ICM, the influence of cooling and tidal effects 
may be needed to explain the details of the observations. 

We discuss the hydrodynamical drag as a heat source for the ICM but conclude 
that it is not likely to play an important role, especially not in stopping cooling flows. 

Key words: galaxies: spiral - galaxies: evolution - galaxies: ISM - galaxies - indi- 
vidual: NGC 4388 - intergalactic medium 



1 INTRODUCTION 



unn fc GottI (j 19721 ) have been the first to propose the 
mechanism of ram pressure stripping (RPS), i.e. the removal 
of a galaxy's gas disc due to its motion through the intra- 
cluster medium (ICM). Since then, this process has been 
studied theoretically and observationally. 

On the theoretical side, several groups have per- 
formed hydrodyna mical simulations, using either SPH 
or grid codes (e.g. 'Abadi et al.l 

Schulz & Struck 2001; Vo Umer et al - -r^-^^^asB,^^^ 

2003; R oediger fc Hensled [20051 ; iRoediger fc Bruggenll2006l ; 



1999; 'Quihs et all l200d; 
2001; Marcolini et all 



Roediger et al.l I2OO6I ). In all these simulations, the model 



galaxy was exposed to a constant ICM wind in order to iso- 
late the ram pressure effect. These studies showed that the 
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ram pressure can remove a significant fraction of a galaxy's 
gas disc or can even strip it completely. In many situations, 
the amount of gas lost cari_be estimate d by a simple ana- 
lytical estimate based on Gunn & Got^j (Il972l ) which com- 
pares the ram pressure to the galaxy's gravitational restoring 
force. 

However, when galaxies move through clusters, they do 
not experience constant ram pressures , but th e ram pressure 
varies along the orbit. IVollmer et all l|200ll ) were the first 
to simulate galaxies in a varying ram pressure. This group 
used a sticky particle code to model the galactic gas disc and 
added the ram pressure as an extra acceleration on all gas 
particles exposed to the win d. This code was applied to sev- 
eral indi vidual gal axies (e.g.lV oUmer e t al.lll999l.l2000l.l200ll: 
Vollmd[20 03; Vo llmer fc Hu chtmcier l2003l : IVollmer et al.l 
2004 120051 . i200tf ) in order to disentangle their ram pres- 
sure histories. In these simulations, the ram pressure was 
allowed to increase and decrease according to a given tem- 
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poral ram pressure profile but did not vary in direction. Be- 
sides this, the sticky particle code cannot model hydrody- 
namical effects such as i nstabilities that have b een s hown 
to play a role . Rece ntly, iRoediger fc BriiggenI (|2007l ) and 
Ijachvm et al.l (|2007l ) have presented hydrodynamical ram 
pressure simulations of gala xies moving through a galaxy 
cluster. Using an SPH code, IJachvm et all (|200'i1 ) have fo- 
cused on rather compact clusters (similar to the Virgo clus- 
ter), so they modelled the ICM-ISM interaction in the inner 
140 kpc of the cluster only. Moreover, in their simulations, 
the galaxy always moves face-on on a strict radial orbit. In 
these compact clusters, the ram pressure peaks are rather 
short and ram pressure stripping becomes a distinct event. 
The ram pressure peak can even be shorter than the strip- 
ping timescale, i.e. the timescale needed to remove the gas 
from the galaxy's poten tial. In such cases, th e galax y loses 
less gas than predicted. iRoediger fc BriiggenI l|2007l ') (here- 
after paper I) have studied RPS in two clusters, one com- 
pact (though not as compact as Virgo) and an extended 
one (similar to the Coma cluster). In this work, the galaxies 
move on different orbits and and also with different incli- 
nations. Here, the ram pressure peaks are not short enough 
that the delay in gas loss becomes important. As long as the 
galaxy does not move close to edge-on or experiences signifi- 
cant continuous stripping, the analyti cal estimate gives good 
predictions of the stripping efhciencv. lTonnesen et al.l (|2007l ') 
have performed a cosmological cluster simulation where they 
can resolve ram pressure stripping with ~ 3 kpc. They find 
that RPS is the main gas loss mechanism in clusters. They 
stress that at a given cluster-centric radius, galaxies can ex- 
perience a variety of ram pressures e.g. due to ambient mo- 
tions in the ICM. From an analysis of cosmological simula- 
tions, |Brfiggen|^|^e^Tuci3 (2007) find that the majority of 
all cluster galaxies experience strong ram pressures within 
their life-times. 

Clearly, RPS provides an enrichment process for the 
ICM, which is known to have a metallicity of about 1/3 solar. 
The metal enrichment o f the I C M has been modelled e .g. by 
ICoral l|2006h : IValdarninil(|2003l ;i: lTornatore et al.l(|2004l ). The 
effect of RPS on the m etal distribution in the I CM has been 
studied numerically bv lDomainko et"al ] (l2006l). which relies 
on the stripping criterion of ICunn fc GottI (|l97l l to model 
the gas loss from the galaxies. Other enric hment processes 
are ga lactic winds ijKapferer et al.l |2006| ). outflows from 
AGN (jMofl et al.ll2007l l. intracluster stars (jZaritski et al.l 
I2OO4I ') . The recent generation of X-ray observatories provides 
information a bout the distribution of metals in the nearby 
clust e rs (e.g. ISchmidt et all I2OO2I: [Matsushita et al.1 l2002l . 



ri 17 



X-ray (jWang et al.1 12004': 'Sun & Vikhlinin"2005'; 'Sun et aL| 
20o4l o r and Ha tails (Gavazzi et al. 2001; Sun ot al. 2007 



Yagi et aI]l2007^ . However, the interpretation of these obser- 



I2OO3I : IChurazov et al.ll2003l : [Pe Grandi et al.ll2004h , which 
the models of metal enrichment have to explain. Recently, 
even the enhanced metallicity in the tail of the ram pressure 
stripped elliptica l galaxy NGC 7619 in the Pegasus group 
ijKim et al.ll2007l ) has been reported. 

Galaxies that experience ram pressure stripping are 
expected to have truncated gas discs but undisturbed 
stellar discs. Such cases have bee n obser v ed: e.g. NGC 
4522 liKeimev fc KoopmannI Il999l . I2001I: iKennev et all 
l200JVollmer et al.ll2004ll. NGC 4548 (jVollmer et al.lll99"i ) 
and NGC 4848 (jVoUmer etHI l200ll ). Deep HI observa- 
tions also revealed long, one-sided tails for se veral galaxies 
(jOosterloo fc van Gorkomll2005l . IChung et al.l [20071 In ad- 
dition to HI tails, there are only very few galaxies that have 



vations as RP S tails is not straigh tforward. E.g. the galaxy 
tails found by IChung et al.l (|2007l ) all belong to Vir go spi- 
rals which are located at projected cluster-centric distances 
of 0.6 to 1 Mpc. Given that ram pressure stripping is ex- 
pected to be strongest near cluster centres, this was a sur- 
prising result. The X-ray and Ha tails tend to be rather 
long (several 10 kpc), narrow (< 10 kpc) and strai ght. More- 
over, especially the X-ray tails seem to be rare: ISun et ah! 
|2003) have searched for additional cases in the Chandra 
and XMM data of 62 galaxy clusters and did not find any. 
Also the HI tails seem to be less common than expected: 
IVollmer fc Huchtmeieil (^007} have carried out deep HI ob- 
servations in the vicinity of 5 RPS candidates in the Virgo 
cluster and did not find more HI than was already known. 

In order to understand the observations of galactic 
wakes, a better theoretical understanding is needed. Here, 
we analyse the simulations presented in paper I, i.e. hydro- 
dynamical simulations of RPS of a disc galaxy along orbits 
in galaxy clusters, with respect to galaxy wakes. We focus 
on the following aspects: 

• Structure of galactic tails: width, length, substructure, 
velocity structure. 

• Comparison to observations 

• Where in the cluster does the galaxy deposit stripped 



gas.'' 



How much heating can RPS provide for the ICM? 



2 METHOD 

Here we analyse the same simulations as in paper I. We 
model the flight of a disc galaxy through a galaxy cluster. 
The galaxy starts at a position ~ 1 to 2 Mpc (depending on 
orbit) from the cluster centre with a given initial velocity. 
We use analytical potentials for the galaxy and the clus- 
ter, as this reduces computational costs significantly. Given 
the high velocities of cluster galaxies, the tidal effect on the 
galaxy is expected to be small. The work o f Moore et al. 
(1996, 1998, 1999). iMastropietro et all (|2005l ) demonstrated 
that only harassment, i.e. the cumulative effect of frequent 
close high velocity encounters between cluster galaxies and 
the overall tidal field of the cluster affects cluster galaxies se- 
riously. Also less massive galaxies are affected more strongly 
by tidal forces. For group environments, there is evidence 
that the tidal forces also influence the ram pressure strip - 
ping efficiency (jMastropietro et al.ll2005l : lMaver et~alll2006l ) . 
We have, however, chosen a massive galaxy (rotation veloc- 
ity 200kms~^) and follow only the first orbit in a cluster 
environment. Thus we expect that tidal forces play a sec- 
ondary role in our case, al though there may b e some in- 
fiuence. Similar to paper I, [jac^metaD presented 
SPH simulations of RPS of a galaxy orbiting through clus- 
ter. They modelled the ICM-ISM interaction only near the 
cluster centre, but included the mutual gravity of all gaseous 
and non-gaseous galactic particles as well as the gravitation 
of the cluster potential. They did not find a significant in- 
fiuence on the stellar and DM components. Thus, we expect 
that our treatment yields reasonable results. 

The orbit of the galaxy is determined by integrating the 
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motion of a point mass through the gravitational potential Table 1. ICM parameters: Core radius, i?icM, beta-parameter, 

of the cluster. In the course of the simulation, the position central ICM density, picMO, and ICM temperature, Ticm, 

of the galaxy potential is shifted along this orbit. for clusters CI, C2 and C 3. Also two parameter sets for the 

Virgo cluster are giv en: from lMatsumoto et al.l l|200Cll ') (MOO) and 

IVollmer et al] l l200j) (VOl) 



2.1 Code 

The sim ulations were performed with the FLASH code 
ijFrvxell et al. 2000) version 2.5, a multidimensional adap- 
tive mesh refinement hydrodynamics code. It solves the Rie- 
mann problem on a Cartesian grid using the Piecewise- 
Parabolic Method (PPM). The simulations presented here 
are performed in 3D. The gas obeys the ideal equation of 
state with an adiabatic index of 7 = 5/3. The size of the 
simulation box is chosen such that the galaxy's orbit during 
the simulation time (3 Gyr) fits into the grid. Depending 
on the orbit, the size of the simulation box ranges between 
(2Mpc)^ and 2 x 5 x 2Mpc^. All boundaries are refiecting. 

The coarsest refinement level has a resolution of Ax = 
62.5 kpc. For most runs, we use 8 levels of refinement, i.e. the 
best resolution is Ax = 0.5 kpc. In addition to the standard 
density and pressure gradient criteria, our user-defined re- 
finement criteria enforce maximal refinement on the galactic 
gas disc and enforce stepwise de-refinement with increasing 
distance to the galaxy. Details are described in paper I. The 
important point for the investigation of the wakes is that we 
limit the refinement outside 50 kpc around galaxy centre to 
Ax = 2 kpc (6 refinement levels), and we limit the refine- 
ment outside 150 kpc around galaxy centre to Ax = 15.6 kpc 
(3 refinement levels). Together with the standard refinement 
criteria, this leads to a typical resolution of Ax = 2 kpc 
at 150 kpc behind the galaxy. Given our resolution restric- 
tion, our analysis and discussion of the wake properties is 
restricted to the first 150 kpc behind the galaxy. A discussion 
of the influence of the resolution on our results is given in 
App.|Xl we flnd our results not to be sensitive to resolution. 

The FLASH code offers the opportunity to advect mass 
scalars along with the density. In order to be able to iden- 
tify the galactic gas after it is stripped from the galaxy, we 
utilise a mass scalar, /, to contain the fraction of galactic 
gas in each cell. Initially, this array has the value 1 in the 
region of the galactic disc and elsewhere. As a result, ev- 
ery cell where / > contains a certain amount of gas that 
has originally been inside the galaxy. At each time-step, the 
quantity fp gives the local density of galactic gas. In the 
context of this paper, we will refer to this gas as galactic gas 
or ISM even if it has left the galaxy. 



2.2 Model galaxy 

The galaxy model is the same as in iRoedieer fc BriiggenI 



-RiCM/kpc 



PiCMo/(gcm 3) Ticm/K 



l|2006r and paper I, i.e. a massive spiral with a flat rotation 
curve at 200 km s"^. It consists of a dark matter halo (1.1 ■ 
10"Mq within 23 kpc), a stellar bulge (10^°Mq), a stellar 
disc (10^^ Mq) and a gaseous disc (5 • IO^Mq). AU non- 
gaseous components just provide the galaxy's potential and 
are not evolved during the simulation. For a description of 
the individual components and a list of parameters please 
refer to RB06. 
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Figure 1. Density and pressure profiles for model clusters CI 
and C3. The profiles for cluster C2 are a factor of 2 lower than 
the ones of cluster CI. For comparison, two profiles for the Virgo 
cluster from the literature are shown. 



2.3 Cluster model 

The cluster is set in hydrostatic equilibrium, where the den- 
sity of the ICM follows a /3-profile, 

-3/213 



picmir) = picMO 



1 + 



RiCM 



(1) 



while the ICM temperature is constant. Given the density 
and temperature distribution and thus pressure distribution 
throughout the cluster, the underlying gravitational accel- 
eration due to the cluster potential is calculated from the 
hydrostatic equation. Here we present simulations for three 
clusters. The parameters are listed in Table [T] along with 
two parameter sets for the Virgo cluster from the literature. 
Figure [T] shows the density and pressure profiles of clusters 
CI and C3. The only difference between clusters CI and C2 
is a factor of 2 in density. These two clusters are compact, 
the ICM density is strongly peaked. Thus, they are simi- 
lar to the Virgo cluster, but not as extreme as Virgo (see 
Fig . [H. Cluster C3 r esembles the Coma cluster (parameters 
see lMohr et al.lll999l ). which is very extended. 



© 2007 RAS, MNRAS OOO.HEs] 



4 E. Roediger and M. Briiggen 



2.4 Galaxy orbit 

The orbit of the galaxy determines its ram pressure history. 
We aimed at constructing orbits with medium to high ram 
pressure peaks, i.e. where the galaxy is expected to lose a 
significant fraction of its gas or is even stripped completely. 
We concentrate on orbits of galaxies that could be regarded 
as falling into the cluster for the first time, i.e. they start 
from a sufficiently large distance from the cluster centre. 
Additionally, the galaxies are bound to the cluster. 

Here we present simulations for four orbits that are sum- 
marised in Fig. [5] Three of these orbits are highly radial, 
whereas the fourth has a large impact parameter. The labels 
of the runs represent the cluster (C-|- number) , small or large 
impact parameter (SM or LG), fast or slow galaxy (FST or 
SLW), inclination (F for near face-on, M for medium, E for 
near edge-on; FE for first near face-on but then near edge- 
on, etc.). In all cases, the galaxy orbits in the j/-z-plane. 
Figure [3] gives a sketch of the orbital plane and galaxy orbit 
in the simulation box. Table [5] lists the simulation runs and 
the amount of lost gas according to paper I. 



3 RESULTS 

As the galaxy moves through the cluster, the ram pres- 
sure stripped gas forms a tail behind the galaxy. Some 
slices in the orbital plane showing the colour-coded gas den- 
sity can be found in paper I and in Fig. IA2I ISM densi- 
ties in the tail are around lO"^'^ gcm""^ close to the galaxy 
(~ 20kpc distance from galaxy centre) and decrease to a 
few 10~^®gcm~^ at larger distances (50kpc). 



3.1 Projected densities 

Figures [S] and [S] show snapshots of the projected ISM 
density for several runs. Each figure is for a different orbit. 
In Figs. 13 and [S] the columns compare cases with different 
inclination but identical orbit. Figure [S] compares two differ- 
ent lines-of-sight (LOS) for the same run. The bottom rows 
of Figs. 13 and [S] summarise the temporal evolution of the 
tail width as a function of distance to galaxy (see discussion 
below). 

The tails of stripped gas stretch along the galaxy's or- 
bit. Three features are prominent in nearly all snapshots: 
the tails show a flaring width, they oscillate along the orbit, 
and are densest close to the galaxy. Details of the structure 
depend on the galaxy's current inclination and the stripping 
stage. In the following subsections, we describe different as- 
pects of the tails. 



3.1.1 Tail length 

Typical column densities are a few 10^^ cm"'^ near the 
galaxy and ~ 2- 10^* cm~^ at large distances (> 100 kpc). A 
typical column density sensitivity limit for current HI obser- 
vations is 10^^ cm~^. Thus, we will define the length of the 
tail as the extent of the light blue contour (~ 2 • 10^^ cm~^) 
behind the galaxy. In this sense, a typical tail length is 
40 kpc. Tails of this length can be found at early times 
(< 1 Gyr) of simulations C1-LG-FST-. . . (two top rows of 
Fig. El), in runs C1-SM-FST-. . . (Fig. P and at early times 
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Figure 2. Summary of galaxy orbits for our simulations. The 
top panel shows the orbits, which are in the y-z-plane. Crosses 
mark the position of the galaxy in intervals of 250 Myr. The next 
panels show the temporal evolution of ram pressure, pram, galaxy 
velocity, Wgai, ICM density along orbit, piCM, distance to cluster 
centre, d, and ratio between ram pressure and local ICM pressure, 
Pram/piCM- For an explanation of the labels see text. 
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Table 2. List of runs. Run labels code the cluster (C+number), small or large impact parameter (SM or LG), fast or slow galaxy (FST 
or SLW), inclination (F for near face-on, M for medium, E for near edge-on; FE for first near face-on but then near edge-on, etc.). Initial 
galaxy coordinates Xg^io and Zg^io are always zero, j/gaio is given in the second column. The third column lists the initial galaxy velocity. 
The fourth and fifth columns list the impact parameter and velocity, respectively. The inclination listed in the sixth column is the angle 
between the galaxy's rotation axis and the j/-axis. The seventh column states the fraction of lost gas according to paper I. 
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Figure 3. Sketch showing the position of the orbital plane in the 
simulation box. Throughout this work, we use two lines-of-sight 
(LOS) to view the galaxy: along the x-axis (towards negative 
direction), which means the LOS is perpendicular to the orbital 
plane; and along the 2-axis (towards negative direction), which 
means the LOS is in the orbital plane. 

(< 1 Gyr) of run C3-SM-FST-MF. We wish to draw at- 
tention to the fact that at these moments the galaxy is 
still far, namely 400 kpc to 1 Mpc, from the cluster cen- 
tre. In the Coma-like cluster C3 (run C3-SM-FST-MF) , at 
t = 500 Myr, the galaxy is even 1800 kpc from the cluster 
centre, and still shows a 50 kpc long tail. Shorter tails are 
found during the first Gyr of runs C2-SM-SLW-. . . (Fig. 
The longest tails are found during peri-centre passage of the 
same run, they extend to beyond 150 kpc. Also the snapshots 
at 1.5 Gyr of runs C1-LG-FST-. . .reveal tails of ~ 100 kpc 
length, although - due to the large impact parameter of this 
orbit - also in this case the galaxy is still more than 400 kpc 
from the cluster centre. 

The length of the tail is influenced mainly by two pa- 



rameters: the ram pressure and the galaxy's velocity. The 
ram pressure determines how much gas the galaxy loses, 
whereas the orbital velocity determines over which volume 
the lost gas is spread. Clearly, high ram pressures lead to 
a high mass loss rate. However, as high ram pressures are 
usually associated with high velocities, the stripped gas is 
also distributed over a larger orbital length or volume. These 
two parameters can be combined into the mass loss per or- 
bital length, which is the crucial quantity that characterises 
many tail properties. Figure [7] displays the evolution of the 
remaining gas disc mass for all runs as a function of dis- 
tance covered by the galaxy. The slope of these functions is 
the mass loss per orbital length. We fltted them piecewise 
with linear functions. Each piece is labelled with its slope 
in lO^A/Q/kpc. The mass loss per orbital length ranges be- 
tween a little below lO'^MQ/kpc to about 7- lO^MQ/kpc. 
Interestingly, this quantity is largest in runs C2-SM-SLW- 
. . . during the peri-centre passage, although the ram pressure 
along most parts of this orbit is smaller than for all other 
orbits. Even during peri-centre passage, the ram pressure for 
this orbit is the second lowest. Here, the low orbital velocity 
causes a high mass loss per orbital length. The case of the 
Coma-like cluster C3 illustrates the opposite extreme of the 
same issue: Here the ram pressure is always higher than in all 
other runs, but the high orbital velocity of the galaxy leads 
to only a medium mass loss per orbital length. As expected, 
we find that the mass loss per orbital length correlates with 
the tail length (see Fig. [S]! . However, additional effects like 
the tail width play a role and lead to some scatter in this 
relation. Also a projection non-perpendicular to the galaxy's 
direction of motion can enhance the projected density of the 
tail while making it appear shorter. 



3.1.2 Tail width 

According to sim ple analytical estimates 

(|Landau fc Lifschit j Il99ll ). the width, w, of the wake 
behind a body moving through a fluid scales with distance 
to the body, d, as y/d in the case of a laminar flow and as 
d^/^ in the case of a turbulent flow. In our simulations, the 
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C1-LG-FST-... C1-SM-FST-... C2-SM-SLW-... C3-SM-FST-... 
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Figure 7. Evolution of gas disc mass as a function of distance covered by the galaxy (top row). Each column is for one orbit (see title of 
column). The colours code runs with different inclinations, see legend. The thin vertical lines mark the peri-centre passage. The bottom 
row displays the relation of time and covered distance. We applied piecewise linear fits to the gas disc mass as a function of covered 
distance. Each piece is labelled with its slope in lO^MQ/kpc. 



flow is clearly turbulent. More precisely, at large distances 
behind the body, the width of the wake scales as 



Fd 



1/3 



(2) 



where F is the drag force the body experiences, p is the 
density of the fluid and U is the flow velocity. The drag 
force can be written as _F" = 0.5cv,ApU^ , where Cw is the 
drag coefficient and A the cross-section of the body. This 
leads to 



50 



kpc(^) 



1/3 



2/3 



lOkpc/ VlOOkpc 



d 



1/3 



(3) 



where we have inserted values typical for our case. Here, r 
is the radius of the remaining gas disc. 

In our simulations, tail widths range between 20 to 50 
kpc at a distance of ~ 25 kpc to the galaxy centre and 30 to 
80 kpc at a distance of ~ 100 kpc to the galaxy centre. Near 
the galaxy, the tail width is similar to the galaxy's cross- 
section with respect to the ICM wind direction. In an order- 
of-magnitude comparison, the simulation results match the 
analytical expectation. However, a direct application of the 
analytical scaling relation to our simulations is not possible, 
as in our simulations many assumptions that went into the 
analytical relation are not fulfilled: The flow past the galaxy 
is not static but, both, flow density and velocity change, the 
body (the galactic gas disc) varies in size. Furthermore, here 
we are interested in the wake close to the body and not far 
away from it. Thus it is not reasonable to fit the d^^'"^ law to 
the galaxy's tail. 

Nonetheless, in order to quantify the flaring ratio of the 
tails, we applied linear fits to the tail widths as a function 
of distance to the galaxy, where examples are given in the 
bottom rows of Figs. |4] and [5] We fixed the j/-axis cut to the 
galaxy's width and fitted the slope to the data. This slope 
is the fiaring ratio. 

For each run, we have derived the fiaring ratio for two 



lines-of-sight: once along the grid's a;-axis and once along 
the grid's z-axis. The a;-axis is always perpendicular to the 
galaxy's orbit, whereas the z-axis is in the galaxy's or- 
bital plane. Figure |5] demonstrates that the flaring ratios 
of the tails are generally independent of the direction of 
the line-of-sight. The only exceptions are some of the early 
{t = 500 Myr) snapshots, where the galaxy's remaining gas 
disc is still large. Figure [TU] displays the dependence of the 
flaring ratio on several quantities. For high ram pressures 
and high galaxy velocities, the flaring ratios are rather small, 
which reflects that for these cases the gas disc is already 
stripped heavily and thus is small. For smaller velocities 
and ram pressures, the fiaring ratios show a large scatter. 
In addition to this, figure [TU] displays the dependence of the 
flaring ratio on the galaxy's inclination with respect to the 
ICM wind direction (top panel) and to the galaxy's current 
diameter (second panel). Both plots reveal correlations in 
the sense that galaxies moving near face-on (small inclina- 
tion angle) and galaxies with large diameters produce tails 
with stronger flaring. The only "exceptions" are again the 
early cases {t = 500 Myr). Combining these two correlations, 
we plotted the flaring ratios over the galaxy's cross-section, 
A — Trr^ig;. cos(j), with respect to the ICM wind direction, 
where rdisc is the radius of the remaining gas disc. This 
plot reveals a clear correlation that small cross-sections lead 
to little flaring and vice versa. A galaxy can have a small 
cross-section either if it is already heavily stripped or if it is 
moving near edge-on. The latter cases are especially inter- 
esting. If such a galaxy is seen perpendicular to the disc, the 
tail is broader than if it is seen edge-on. However, the flar- 
ing for both line-of-sights is small. This reflects the known 
dynamics of turbulent wakes: according to Eq. [S] at a given 
distance, d, behind an object, the flaring of the wake scales 
as ^ oc A^^^ , i.e. the wakes of larger objects show stronger 
flaring. As emphasized above, however, several assumptions 
that went into Eq.[3]are not fulfilled in our simulations, thus 
there is little hope of recovering the A^^^ relation. We only 
flnd that a larger galactic crossection leads to stronger flar- 
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Figure 4. Projected galactic gas density for runs Cl-LG-FST- 
EF (left column) and Cl-LG-FST-FE (right column; same or- 
bit, different inclination). Lowest contour is 10~®g cm~^=6.0 ■ 
10^^ cm"'^. The contour spacing is half an order of magnitude. 
The time of each panel is denoted in its upper left corner. The 
label "axis: x" means the projection is done along x-axis, which 
means perpendicular to the galaxy's orbital plane in our simula- 
tions (see Fig.O. The galactic orbit is marked by the black line. 
The plus-sign marks the galactic centre. The coordinates at the 
axes are in the cluster-centric system. The bottom row displays 
the tail width as a function of distance to the galaxy at different 
time-steps as measured in the snapshots. Same colours are for the 
same time (see legend), thin lines are for projection along grid's 
0-axis, thick lines for projection along x-axis (like this Figure). 



ing. Nonetheless, these characteristics suggest that the flar- 
ing is determined by the turbulence in the ICM flow past 
the galaxy. 

Effects commonly seen in simulations of harassment 
and tidal interactions arc tilts an d warps of the galac - 
tic discs (e.g. Mooro ct al. 1998; Ma stropietro et ahllioosh . 
These may change the orientation of the galactic disc, i.e. 
stars and gas, continuously. This differs from what we have 
assumed here. However, the effect of the cluster's tidal field 
is strong only near the cluster centre, where the gas disc 
is already heavily ram pressure stripped. Of course interac- 
tions with other cluster galaxies may modify the gravita- 
tional potential of both galaxies and thus also influence the 
ram pressure stripping efficiency and characteristics of the 
wake. In this work here, we have isolated the effect of ram 
pressure stripping. 

The description of the tail widths above considered the 
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full tails. If we again adopt a column density limit of ~ 
10^^ cm~^ (light blue contour), the tail flaring is much less 
pronounced, often not even detectable. Generally, then the 
tail widths are similar to the cross-section of the remaining 
as disc. 



3.1.3 Mass distribution along the tail 

Fig. [TT] shows the galactic gas mass per orbital length as a 
function of distance to the galaxy, D, for three exemplary 
cases: we have integrated the projected gas densities perpen- 
dicular to the galaxy's current direction of motion. As the 
orbits are only slightly curved, for the 150 kpc behind the 
galaxy, this corresponds to very good approximation to an 
integration perpendicular to the orbit. The peak at D — 
is the gas still trapped in the galaxy. With increasing D, 
the mass per orbital length decreases and often saturates 
at f ~ 70 kpc. However, this general behaviour is superim- 
posed with substantial substructure. Especially in the tails 
at early simulation times (t = 500 Myr), the mass per orbital 
length does not saturate but continues to decrease. 

The gas mass per orbital length should record the gas 
loss history of the galaxy. In Fig. Owe have applied piecewise 
linear fits to the evolution of the galaxies' gas disc masses as 
functions of covered distance. The slopes give the mass loss 
per orbital length and thus should be the same as the gas 
mass per orbital length in the tails. To compare both quanti- 
ties, we have plotted the mass losses per orbital length as de- 
rived from Fig.[7]as thin horizontal lines in Fig. 1111 where the 
colours code the time-steps according to the legend. Indeed, 
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Figure 11. Galactic gas mass per orbital length as a function of distance to the galaxy for three exemplary runs. Each panel is for 
one run (see legends), the line colours code different time-steps, see also legends. The thin horizontal lines mark the predicted mass per 
orbital length as derived in Fig. [7] 



when the mass per orbital length saturates, it saturates at 
approximately the level marked by the thin horizontal lines. 
There are several reasons for the differences between the de- 
rived mass loss per orbital length and observed gas mass per 
orbital length; First of all, the mass loss per orbital length is 
no constant but also varies on short time- and length-scales. 
Secondly, the stripped gas has to be accelerated away from 
the galaxy, which also influences the gas distribution along 
the tail. Moreover, gas can be temporarily trapped in tur- 
bulent eddies, which places it in a position that cannot be 
predicted by the mass loss and acceleration. 

3.2 Velocities in the wakes 

Some slices in the orbital plane showing the gas density and 
the velocity field can be found in Fig. \Mi The ICM is flow- 
ing around the gas disc. In cases where the galaxy moves 
supersonically, a bow shock is prominent. In the wake, the 
velocities are generally smaller than the ICM wind velocity. 
Additionally, the velocity field shows a turbulent structure. 

Figures [T5] to [TJ] summarise the information about the 
velocity in the wakes for some exemplary cases. For an ISM- 
density-weighted random subset of grid cells we calculate 
Wpar, the velocity component anti-parallel to the galaxy's 
direction of motion, («i:gai, "Ut/gai, «zgai). We also calculate 
iiperp, the component perpendicular to Vpar. By definition, 
'^pcrp is always positive. Both velocity components are plot- 
ted as a function of distance to the galaxy, D, i.e. the projec- 
tion of each grid cell's position vector in the galactocentric 
frame onto (— w^gai, — ^^ygah — I'zgai)- Additionally, we show 
the velocity component along the grid's s-axis, Vx, which in 
our simulations is always perpendicular to the galaxy's or- 
bital plane. All velocity components are given in the galaxy's 
rest frame, where the galaxy's motion translates into an ICM 
wind fiowing past the galaxy. 

3.2.1 Velocity component anti-parallel to iTgai 

The overall structure of the Upar-plots is the same for all 
runs and all time-steps. The gas in the galactic disc appears 
as a blob between D = ±20kpc and Upar ~ ±200kms~^, 
which is due to the disc's rotation. With increasing distance 
to the galaxy, Upar increases towards the ICM wind velocity, 
which refiects the acceleration of the stripped gas by the 
ICM wind. At a distance of about 80 to 100 kpc behind the 



galaxy, the acceleration of the stripped gas is completed, 
Vpa,i does not increase any further. For a fixed distance to 
the galaxy, the velocity range in the tail is large. Inside the 
galactic disc, a width of 400kms~^ is expected due to the 
rotation. This behaviour is also refiected in the Wpar-plots. In 
the region where the stripped gas is accelerated, the velocity 
width is larger than this, it can reach 1000 km s~^. Beyond 
distances of ~ 40 kpc behind the galaxy, the velocity width 
decreases. At distances beyond 100 kpc, the velocity width 
is again 400 to 500kms~^ and remains constant. 

We have marked the current ICM wind velocity (or or- 
bital velocity of the galaxy) in each plot by a horizontal line. 
Only very few grid cells reach Upar above this line, which is 
not surprising, as the ICM wind cannot accelerate gas be- 
yond its own velocity. The acceleration of the stripped gas 
always proceeds in a rather similar fashion. We have plotted 
the same dashed diagonal line in all plots in order to aid the 
eye to compare the slope of iipar between different snapshots. 
The slope of the dashed line is representative for snapshots 
with ram pressures around 10~^^ ergcm"^. For stronger ram 
pressures, which are also associated with higher orbital ve- 
locities, the slope is somewhat larger, for weaker ram pres- 
sures somewhat smaller. Especially for the later time-steps 
of run C3-SM-FST-MF (Fig.[Tl|, where the orbital velocity 
is much larger than in the other cases, the velocity slope 
along the tail is steeper. The large velocity width in the tail 
near the galaxy is caused by a superposition of several pro- 
cesses. Firstly, also here the velocities originating from the 
rotation of the gas disc add to the velocity width. Secondly, 
denser ISM clouds need a longer time to be accelerated, 
while less dense structures are accelerated more easily. A 
third effect is the turbulence in the wake. As in our simula- 
tions the spatial resolution beyond 150 kpc from the galactic 
centre is limited to 15 kpc, these simulations cannot provide 
reliable information beyond this distance. 

We observe some occasional backfall of stripped gas 
well before the ram pressure peak. The backfall is tem- 
poral and becomes evident as negative Wpar closely behind 
galaxy (e.g. early panels of Fig. I13p . and as local max- 
ima/plateaus in temporal evolution of gas mass in disc 
region (see paper I, or Fig. [7]). A corresponding feature 
is evident in earlier simulations with constant ICM wind. 
There, the gas disc mass as a function of time shows a 
local maximum or plateau immediat ely after the instanta- 
neous stripping phase (e.g. grid codes: [Marcolini et al.ll2003l : 



© 2007 RAS, MNRAS 000,[TH23] 



10 E. Roediger and M. Briiggen 



m 
E 



1500 
1000 
500 



1^88 

■fOOO 
500 


m% 

500 

b 
-500 
-1000 



0.4 Gyr 



0,8 (Syr 



1,2 Gyr 



1.6 Gy r 





SIM 



H =H 



H \ 



50 100 1500 50 100 1S0O 50 100 

true dsslarrce to galaxy in Vpz 



1500 



50 100 150 



Figure 12. Three different velocity components in the galaxy's tail, as a function of distance to the galaxy. For run Cl-LG-FST-EF. 
The top row shows the velocity component along the tail. The middle row shows the component perpendicular to the galaxy's direction 
of motion, and the bottom row the component along the grid's x-axis, i.e. perpendicular to the galaxy's orbital plane. Each column is 
for one time-step, as indicated. The velocities are given in the galactic rest frame. An ISM density-weighted random subset of grid cells 
is shown. 

The solid horizontal lines in the top row panels mark the current ICM wind velocity, i.e. the negative of the galaxy's current velocity. 
An identical diagonal dashed line is plotted in the top row panels in order to aid the eye to compare the slope of the velocity gradient. 
The horizontal lines in the middle row panels mark and 200kms~^, the horizontal lines in the bottom row panels mark ±200kms~^. 




iRoediger fc Henslerl l2005l: iRoediger fc BruggenI l2006l . but 
also in SPH: ISchulz fc Struckll200ll '). Of couse. this feature 
is only observable if the evolution of the model galaxy was 
followed for a sufficiently long time. However, this backfall is 
not observed in all runs, but mainly while the galaxy moves 
near face-on. While the galaxy moves near edge-on, hardly 
any backfall is observed. Thus we conclude that the tempo- 
ral backfall in our simulations is caused by the turbulent flow 
around the remaining gas disc. Also the shape of the galaxy's 
total potential supports such this backfall behaviour, as, in 
the face-on case, the ram pressure meets the strongest grav- 
itational restoring force a few kpc behind the disc plane 
(see e.g. ISchulz fc Struck I2001I : IRoediger fc Hensled l2005l : 



Ijachvm et al ] I2OO7I . but it is also the case for our model 
galaxy). 

3.2.2 Velocity component perpendicular to Ugai 

Velocities perpendicular to galaxy's direction of motion, 
Wpcrp, can be due to the disc's rotation as well as the ICM 
flow around the galaxy. Moreover, turbulence in the wake 
adds to Uporp. The component Wporp causes the tail flaring. 

If all three sources were responsible for the tail flaring, 
we should observe that in cases where the remaining gas 
disc is still large but the galaxy is moving edge-on, the flar- 
ing differs between LOS perpendicular and parallel to the 
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Figure 14. Same as Figs. [121 but for run C3-SM-FST-MF. 



galactic disc. If such a galaxy is seen face-on, the rotation 
and the ICM flow around the galaxy contribute to the veloc- 
ity component perpendicular to the LOS and tail direction. 
In contrast, if such a galaxy is seen edge-on, neither ICM 
flow nor rotation contribute to the velocity component per- 
pendicular to the LOS and tail direction. Thus, one could 
expect that the tail flaring should appear stronger in the 
case where the galaxy is seen face-on than in the case where 
the galaxy is seen edge-on. However, we observe similar flar- 
ing ratios for both LOS. Thus we can conclude that mainly 
turbulence causes the flaring. 

The velocities perpendicular to the galaxy's direction of 
motion range between and up to 700kms~^. In near edge- 
on cases, Vperp ranges only between and 500kms~^. The 
range of Vperp is also smaller for small ram pressures and 
galaxy velocities. This behaviour corresponds to the rate 
of tail flaring. Again we observe that in cases where the 
galaxy's cross-section is large, the wake is more turbulent. 

The velocity component along the grid's a;-axis, Vx, is 
always perpendicular to the orbital plane and thus perpen- 
dicular to the galaxy's direction of motion. However, un- 
like Uperp, it does not only contain information about the 
amplitude of the perpendicular velocity components, but 
also directional information. It often shows an oscillating 
behaviour, which causes the oscillation of the galactic tails 
along the orbits as described in Sect. 13. ll These oscillations 
resemble von-Karman vortex streets. 

The amplitude of the turbulence in the wakes seen in 
our simulations is comparable to the ICM turbulence gen - 
erated by o ther proces ses: e.g. Norman fc BrvanI l)l999l l. 
iDolag et all f2005 l and IVazza et all (|2006l l find turbulent 
velocities due to structure formation and subclump infall of 



some 100kms~ MKiml lj2007h find that gravitational wakes 
of galaxies produce turbulence at a level of < 220kms~^. 



3.3 Distribution of stripped material throughout 
the cluster 

The galaxy's gas loss history is reflected in the distribution of 
the stripped gas throughout the galaxy cluster. In Fig. llSlwe 
show this distribution for the orbit C1-LG-FST-. . . for two 
different galaxy inclinations. Clearly, different inclinations 
lead to different gas loss histories and thus differing distribu- 
tions of the stripping gas throughout the cluster. In order to 
describe the distribution of stripped gas along the galaxy's 
orbit, we have calculated the ISM mass inside a sphere of ra- 
dius 150 kpc around each orbit point. Dividing this mass by 
the orbit length inside the sphere (which is 300 kpc) yields 
the local ISM mass per orbital length, averaged over 300 kpc. 
This quantity is plotted in Fig. [TS] for different time-steps 
of simulation run Cl-LG-FST-MF. Figure [IT] displays the 
same quantity for all runs at a fixed time-step t — 2 Gyr. 
We had to use a rather large averaging length due to the 
large wake width at large distances from the galaxy. Using 
a smaller smoothing length would have meant to miss some 
gas. 

Figure [17] again reflects the fact that the amount of gas 
that the galaxy deposits locally depends on its mass loss per 
orbital length. E.g. in the simulation run corresponding to 
the bottom panel, the ram pressure is always larger than in 
the other simulations. However, here the galaxy is also mov- 
ing faster than in all other simulations, and thus spreads its 
gas over a larger volume and larger orbital length. Conse- 
quently, the local amount of stripped gas is not higher than 
in the other simulations. 
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Figure 15. Slices in the orbital plane showing the ISM distribution throughout the cluster colour-coded for different time-steps. The two 
top rows are for run Cf-LG-FST-FE, the two bottom rows for run Cl-LG-FST-EF. The Ihs column displays the density of (stripped) 
ISM, piSM, the rhs column displays the local ISM fraction, Pism/picm- The black line marks the galaxy's orbit. 
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Figure 16. Distribution of stripped gas along orbit at different 
time-steps for run Cl-LG-FST-MF, averaged over 300 kpc. The 
thin black line is the prediction based on the analytical estimate 
of the stripping mass (see paper I). The gray line is the prediction 
based on the numerical result of the bound gas mass as a function 
of covered distance - however — shifted by 150 kpc along the orbit. 
The zero-point of the a;-axis is shifted to peri-centre passage. 



The amount of stripped gas found per orbital length has 
to be closely linked to the galaxy's mass loss rate per orbital 
length discussed in connection with Fig. [T] In our simula- 
tions, we know the amount of ISM bound to the galaxy po- 
tential. If we assume that the stripped gas is deposited into 
the ICM exactly where it is lost, we can predict the expected 
ISM mass per orbital length. The gray line in Fig. 1161 shows 
this function, however, shifted by 150 kpc along the orbit. 
The thin lines in Fig. [17] show the same for the other runs. 
This means that the stripped gas is deposited into the ICM 
about 150 kpc from the position where it was lost from the 
galaxy's potential. In other words, the stripped gas follows 
the galaxy for about 150 kpc before finding its final position 
in the cluster centre. Thus, in our simulations, the stripped 
gas is deposited into the ICM rather locally. We have also 
plotted the prediction for the distribution of the stripped gas 
based on the analytical estimate of the galaxy's stripping ra- 
dius (Gunn& Gott criterion, see paper I). If the galaxy lost 
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Figure 17. Distribution of stripped gas along orbit at t = 2 Gyr, 
averaged over 300 kpc. Each panel is for one orbit, different in- 
clinations are colour-coded (see legend). For each run, the thick 
line is the distribution of stripped gas as seen in the simulation. 
The thin line of the same colour is the prediction based on the 
numerical result of the bound gas mass as a function of covered 
distance — however - shifted by 150 kpc along the orbit. The black 
line is the prediction from the analytical estimate. The zero-point 
of the a;-axis is shifted to peri-centre passage. 



its gas according to this criterion and if the stripped gas 
remained where it was lost, its distribution along the orbit 
would be as shown by the thin black lines in Figs. [16] and 
1171 In most cases, the simulations and the analytical predic- 
tion differ, according to the simulations the stripped gas is 
deposited later along the orbit. 

We need to point out that results regarding the distri- 
bution of the stripped gas throughout the cluster are approx- 
imations, as in our simulations we have limited the spatial 
resolution outside 150 kpc from the galaxy. However, accord- 
ing to the velocity plots, the deceleration of the stripped gas 
is finished within 150 kpc from the galaxy. Moreover, if the 
spatial resolution is improved by a factor of 2 throughout the 
simulation box, the results are very similar (see Figs. I All and 
IA5|) . Thus, we conclude that our results do not suffer from 
insufhcient resolution. 

Typical ISM fractions in the wake ~ 100 kpc behind 
the galaxy are around 20%. Assuming an ISM metallicity 
of about solar and an ICM metallicity of about 0.4 solar, 
this leads to a metallicity of 0.52 in the wake. For an ISM 
metallicity of twice solar and ICM metallicity of 0.2 solar, 
the resulting metallicity in the wake is 0.56 solar. In how 
far the wake is detectable in X-ray metal maps depends also 
on the temperature in the wake. As we do not know the 
temperatures in the wake, we cannot make predictions at 
this point. 

Obviously, RPS is a source of metals for the ICM. The 
cumulative effect of the whole galaxy population on the 



ICM is studied by IPomainko et al. (|2006|) and compare d 
to other enrichment processes by ISchindler et all (|2005h : 
iKapferer et"al] (|200'it : iMoll et all (|200'if ). Here, the enrich- 
ment by RPS is modelled by applying the classical Gunn 
& Gott criterion to the galaxy particles and thus calculat- 
ing a mass loss rate for each galaxy. The gas predicted to 
be stripped is added locally to the ICM. In this treatment, 
galaxies lose gas only as long as the ram pressure is increas- 
ing, only on their way towards the cluster centre. Our sim- 
ulations show that galaxies lose their gas somewhat more 
slowly than predicted, stripped gas is also found along the 
orbit after peri-centre passage. These characteristics could 
lead to a slightly broader distribution of ram pressure in- 
duced metals in the ICM. 



3.4 Energy input into the ICM 

The galaxy carries a large amount of kinetic energy. The 
ram pressure working on the gas disc is associated with a 
drag force that decelerates the galaxy. This mechanism pro- 
vides a heating source for the ICM. The energy lost by the 
galaxy can be thermalised in the ICM either via decaying 
turbulence or via viscosity. Here we want to investigate the 
relevance of this heating mechanism. 

A body of mass, M , and cross-section. A, that is subject 
to the ram pressure, Pram, experiences the drag force, -Fb = 
— J3ramj4 aud the associated deceleration 



Pram ^ 

The energy loss rate for this body is 
d 



At 



-Ekin = —Mva-D = —Upram^. 



(4) 



(5) 



The more relevant quantity in our context is the energy loss 
per orbital length. 



ds 



(6) 



We can compare this quantity to the local turbulent energy 
per orbital length in the galaxy's wake, which is 



F 

ds 



- PiCM IJturb ^wako , 



(7) 



where uturb is the typical turbulent velocity in the wake and 
^wakc the wake's cross-section. Assuming that the loss of 
kinetic energy is converted locally into turbulence, energy 
conservation requires that these two quantities are the same. 
Making use of the relation pram = ^CwPiCMf|aii this leads 
to the following relation between the body's (the galaxy's) 
orbital velocity, Ugai, and the ratio of the cross-sections of 
the body and the wake: 



fturb 
Wgal 



A 



(8) 



As usual, Cw is the drag coefficient that parametrises the 
response of the body in question to a flow, it depends on 
the body's shape and surface properties. For a disc-like body 
we have Cw ~ 1 if the disc is moving face-on. In case of our 
galaxy, the ratio between the gas disc's cross-section and the 
wake's cross-section is ~ 4 to 9. For a typical orbital velocity 
of 1000 km s~^ this leads to turbulent velocities of the order 
of 400kms~^, which is what we observe in our simulations. 
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Now we want to address the question of the relevance 
of this process as a heating mechanism for the ICM. Given 
that the analytical estimate for the stripping radius gives a 
reasonable result for a large range of situations (see paper I), 
we can use A = Trrlij.^, as an estimate for the galaxy's cross- 
section, where raise is the stripping radius derived from the 
analytical estimate. Furthermore, we assume that the kinetic 
energy lost by the galaxy is available for ICM heating locally. 
Thus, we can calculate the local heat gain per orbital length 
to be 



-C/hcat 

as 



as 



(9) 



However, not all of this energy can really be used to heat the 
ICM. Some part is consumed in removing the galaxy's gas 
disc from the galactic potential. The total binding energy of 
the gas disc is ~ 1.5- 10^* erg, and it scales approximately lin- 



ear with disc radius with a slope of 



■ lO'^" erg/ kpc. 



Thus, the amount of energy consumed in this process per 
orbital length is 



ds "^"^ drdiac ds 



8-10 



,56 erg drdis, 
kpc ds 



(10) 



where the change of disc radius per orbital length, 
also given by the analytical estimate. Figure[l8]compares the 
energy input into the ICM due to the drag deceleration of 
the galaxy and the energy lost by the removal of the stripped 
gas from the galaxy's potential for the four different orbits. 
The energy budget is given locally along the galaxy's orbit, 
as energy per orbital length. Clearly, the input of energy into 
the ICM due to the galaxy's drag deceleration dominates. 
The radiative loss from the ICM by thermal bremsstrahlung 
was not considered here as due to its strong dependence on 
ICM density it can be relevant only close to cluster centres. 

In order to judge the importance this energy gain for 
the ICM, we can compare the original amount of thermal 
energy along the wake, eAwako, with this heat gain. The 
thermal energy density, e, can be expressed in terms of the 
sound speed, c^: 



' PicmCs , 



(11) 



where we have made use of the fact that in our case (7 — 
1)7 ~ 1. Thus, the ratio of heat gain and thermal energy is 



d-Eheat /ds _ fgal^dis 
— U.OCvi 



(12) 



Usually, the wake diameter is at least a factor of 2 larger 
than the diameter of the gas disc. In our simulations, the 
orbital velocity reaches Mach numbers of 2 only in the clus- 
ter centre. Thus, the energy gain by ram pressure decelera- 
tion would correspond to a local energy (and temperature) 
increase by a factor of 2 at maximum. In principle, the heat- 
ing could be strongest near the cluster centre as there the 
galaxy moves fastest, but due to ram pressure stripping also 
here the gas disc is rather small or even gone completely. 
Thus, the process of drag deceleration is unlikely to be able 
to stop cooling flows. In summary, this means that ICM 
heating due to ram pressure deceleration of galaxies is not 
expected to play a crucial role. This conclusion is not sur- 
prising if we remember that the ICM temperature as well as 
the galaxies' velocity dispersion are determined by the same 
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Figure 18. Energy gain for the ICM along the orbit. For four 
difTeront orbits (see label in upper right corners). The zero-point 
of the X-axis is shifted to peri-centre passage. 



gravitational potential and the total mass of ICM is larger 
or comparable to the total mass of all galaxies. This means 
that the total amount of thermal energy in the ICM is com- 
parable or larger than the total kinetic energy of all cluster 
galaxies. Moreover, the galaxies lose only a fraction of their 
kinetic energy by drag deceleration. In our calculation we 
have also assumed that the galaxy is moving (near) face-on. 
For galaxies moving near edge-on, the drag force as well as 
the heating rate will be somewhat smaller. 

Another process to convert the kinetic energy of the 
galaxies into thermal ene rgy of the ICM is dynamical fric- 
tion (e.g. iKim et al.|[2005l and references therein). The decel- 
erating force working on a galaxy of mass, Mgai as it moves 
through a homogeneous medium of density, p, with velocity, 

Wgal, is 



47V pG^M^^i 



(13) 



"gal 



According to lOstrikej l|l999l l. the efficiency factor / sum- 
marises the dependence on Mach number and impact param- 
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eters, it is in the range of 1 to 10. The density, p, of the sur- 
rounding medium is the ICM density plus the DM density. 
As an approximation, we will use I — 1 and p — lOpicM, and 
a total galaxy mass of 5-10" Mq. Then we can calculate the 
dynamical friction force, which is equivalent to the galaxy's 
energy loss per orbital length. The result is also shown in 
Fig. 1181 The energy available from dynamical friction is com- 
parable to the energy available from hydrodynamical drag. 
However, dynamical friction is stronger in cluster centres. At 
first glance, also the hydrodynamical drag should be strong 
in cluster centres as the ram pressure is highest there, but 
the high ram pressure causes ram pressure stripping and 
thus small galactic gas discs and, consequently, small drag 
forces. 

Given the deceleration by dynamical friction, aoF = 
Fdf /Afgai , we can estimate how much the galaxy gets decel- 
erated during one orbit: 



of 



100 kpc. If we adopt a column density limit of 



At; 



At 



18kms" 



10-26 gcm-3 



gal 



5- 10" Mq 



.1000 kms- 



At 



(14) 



which implies that deceleration due to dynamical friction is 
negligible as long as only the first orbit is concerned. As the 
strength of the hydrodynamical drag force is similar to the 
dynamical friction, the same applies to the hydrodynamical 
drag. 



4 DISCUSSION 

4.1 Summary of tail properties 

The overall picture of ram pressure induced galactic tails 
as derived from our simulations can be summarised as fol- 
lows (we use the galactic rest frame). There are two main 
processes: 

• The ICM wind accelerates the stripped gas along the 
orbit away from the galaxy. The acceleration is finished ~ 
100 kpc behind the galaxy. 

• Additionally, mainly turbulence in the wake leads to 
a flaring of the tail. The turbulence and hence flaring ra- 
tios depend on the galaxy's cross-section with respect to the 
ICM wind direction: Small cross-sections lead to little flar- 
ing, large cross-sections lead to stronger flaring. 

This leads to the following characteristics: 

• The tails of stripped gas stretch along the galaxy's or- 
bit. They oscillate slightly along the orbit and show a flaring 
width. 

• Local ISM densities in the tail are around 10~^^ gcm""^ 
close to the galaxy (~ 20 kpc distance from galaxy centre) 
and a few 10~^* gcm""^ at larger distances (50 kpc). 

• Projected ISM densities are a few 10^^ cm"^ near the 
galaxy and ~ 2-10^* cm"^ at large distances (50 to 100 kpc). 

• The density and length of the tail (for a given column 
density limit) are set by the galaxy's mass loss per orbital 
length. 

• Widths range between 20 to 50 kpc at a distance of 
~ 25 kpc to the galaxy centre and 30 to 80 kpc at a distance 



10 cm~ , the tail width is similar to the galaxy's apparent 
cross-section with respect to ICM wind direction and does 
not show flaring. 

• For the same column density limit, typical tail lengths 
are 40 kpc. Even in compact clusters like Virgo, such tails 
can also be found for galaxies that are still 0.5 to 1 Mpc 
from the cluster centre. In extended clusters like Coma, 40 
kpc long tails can be found even at cluster-centric distances 
of 1800 kpc. For high mass losses per orbital length (e.g. a 
passage close to the cluster centre with moderate velocity), 
tail lengths can reach 150 kpc. 

• At a distance of ~ 70 kpc, the galactic gas mass per or- 
bital length along the tail converges roughly to the galaxy's 
mass loss per orbital length. Typical values are around 
3 • 10^ Mq/ kpc. The mass per orbital length along the tail 
can also be derived for observed galaxies. In order to disen- 
tangle a galaxy's mass loss history from this quantity, the 
tail needs to be observed to at least this distance behind the 
galaxy. With typical tail widths of ~ 50 kpc, a sensitivity 
limit of a few 10^* cm"^ is required. 



4.2 Limits of this work 

In our simulations, we neglected cooling and thermal con- 
duction. Thus, we are unable to predict the temperature 
in the wakes and the spectral range in which they are ob- 
servable. Rece nt observations of long galactic tails in the 
Virgo cluster (lOosterloo fc van Gorkomllioosl . IChung et al.l 
l2007h suggest that the stripped gas could remain cool for 
a few 100 Myr. In contrast, the long X-ray tail observed 
bv lSun et all l|2006l ) - if it is ram pressure induced - sug- 
gests that at least some part of the stripped gas is heated 
to higher temp eratures. Cooling may a lso affect the strip- 
ping efficiency (|Mastropietro et al.ll2005h . However, in order 
to include cooling in a propper way, also star formation and 
supernova feedback would have to be modelled, i.e. the mul- 
tiphase nature of the ISM would have to be taken into ac- 
count. This is beyond the scope of this paper but will be a 
highly relevant subject of future work. 

Cooling may also influence the morphology of the tails. 
Dense regions may cool and contract further, thus forming 
a clumpy tail. However, for the densities (< lO^^^gcm"'') 
and temperatures (~ 10^ K) in the wakes in our current 
simulations the cooling time is at least several 10 Myr and 
thus comparab le or longer than th e dynamical timescale. 

Moreover, IVietri et al.l l|l997l l have shown that in the 
case of dense cool clouds in moving though the warm ISM, 
radiative cooling can prevent the disruption of these clouds 
by the KH instability. The instability of shear flows at the 
boundary of an object that moves through some fluid are the 
main reason for the generation of a turbulent wake. Thus, 
one may wonder if the presence of cooling could also pre- 
vent the generation of turbulence in the wake of our model 
galaxies. One major difference to the work of IVietri et al.l 
(1997) is, however, the density and t emperature range u nder 
consideration. In the simulations of IVietri et al.l (|l997h . the 
cooling times inside the clouds are always much shorter that 
the dynamical timescale. This is not the case in the ICM, 
here cooling times are at least comparable to the dynamical 
timescale . Moreover, a close inspection of the snapshots pre- 
sented in IVietri et al.l (|l997l ) reveals that the flow of the hot 
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gas still shows some eddies, even when coolin g; is included. 
In a similar context, iVieser fc Hensleil |20o3) have studied 
the influence of thermal conduction on the survival of cool, 
dense clouds. They find that the interplay of cooling and 
thermal conduction can stabilise the clouds against KH in- 
stability. Again, the densities and temperatures studied are 
different from our case. Thus, the influence of cooling and 
thermal conduction on the dynamics of the wakes is another 
important task for future studies. 

We also neglected viscosity in our simulations. The am- 
plitu de of viscosity in th e ICM is still a matter of debate 
fe.g. lRevnolds et alllioosl : iRuszkowski et al.ll2004 ). We note 
that the same applies to the heat conduction. A higher vis- 
cosity prevents turbulent motions in the wakes and thus 
would suppress tail flaring. Consequently, in a viscous ICM, 
ram pressure induced tails would be narrower than the 
ones in our simulations. Moreover, the velocity width across 
the tail would be smaller. A comparison between viscous 
and non-viscous RPS simulations and detailed observations 
could help to measure the viscosity of the ICM. 

Despite the many possibilites to improve our model, 
this model provides an important starting point to compare 
theory and observations. 



4.3 Comparison with other simulations 

4.3.1 Constant ICM wind simulations 

iRoediger et al.l (j2006l ) have studied galactic wakes in ram 
pressure simulations using a constant ICM wind. Also these 
simulations showed the acceleration of the stripped gas away 
from the galaxy and the flaring of the tails. However, in these 
simulations, the tails were generally broader than the ones 
presented here. In static wind simulations, the increase of 
Vpa,T due to acceleration takes somewhat longer. However, 
the constant wind simulations suffered from the difficulty 
that the flow had to be initialised at full strength, which 
is artificial. Concerning the dynamical aspect, the cluster 
crossing simulations presented here are by far more realistic. 



4.3.2 Sticky-particle simulations 

Also the sticky-particle simulations of VoUmer et al. (see 
Sect. [T|) describe RPS of galaxies on cluster orbits. Accord- 
ing to their simulations, regarding the particle density, the 
stripped gas often forms a dense arm that originates at one 
edge of the galaxy and is much narrower than the galaxy's 
cross-section. The low particle density extent of the tail has 
the same width as the galaxy and shows some flaring de- 
pending on stripping geometry. In their work, they usually 
concentrated on the gas distribution close to the galaxy, thus 
there are no predicted HI maps for long tails. 

For galaxies moving near edge-on, the sticky particle 
simulations predict a characteristic backfall of stripped gas 
some time after peri-centre passage. We do not observe such 
a backfall in any of our simulations, which is most likely 
due to the fact that our model clusters and thus ram pres- 
sure peaks are somewhat broader than the cases studied by 
IVollmer et all (|200ll 'l. However, we observe a temporal back- 
fall of gas for near face-on cases while the galaxy approaches 
the cluster centre. This backfall is caused by the turbulent 



velocity structure in the wake and can thus not be present 
in the sticky-particle simulations. 

4.3.3 SPH simulations 

Ijachvm et al.l l|2007h presented SPH simulations of RPS of 
disc galaxies on cluster orbits. Also in their work, internal 
processes in the ISM (cooling, star formation, feedback) were 
neglected; the ISM was treated isothermally, while the ICM 
was treated adiabatically. Additionally, they restricted the 
ram pressure interaction to the inner 300 kpc (in diameter) 
of the cluster. Most of their model clusters were more com- 
pact than our model clusters. This fact and the restriction 
to the inner cluster part may be the reasons why also this 
group observes reaccretion or backfall of stripped gas af- 
ter peri-centre passage, while our simulations do not show 
such a behaviour. The galactic wakes in these simulations 
also seem to differ from the wakes in our simulations. How- 
ever, a detailed comparison is difficult as they mostly show 
only the first 20 kpc behind the galaxy. In the simulations of 
IJachvm et all (|2007l ), the opening angle of the tail is influ- 
enced mainly by the ram pressure strength: when the ram 
pressure is still low and the galaxy moving slowly, the tails 
appear to be flaring strongly, and with increasing ram pres- 
sure and galactic velocity the opening angle decreases and 
remains moderate after peri-centre passage. Apparently, the 
tail width at a distance of 100 kpc behind the galaxy is 
less than 50 kpc, which is not much more than the origi- 
nal diameter of the galaxy. Additionally, the velocity in the 
wake shows mainly a slow motion in ICM wind direction 
and only very small velocity components perpendicular to 
the galaxy's direction of motion. Thus, the tail opening an- 
gle in these simulations seems not to be caused by a sig- 
nificant perpendicular velocity component as in our case. It 
mainly reflects that flrst the outer parts of the gas disc are 
stripped, and, as the galaxy approaches the cluster centre, 
gas from smaller and smaller radii is stripped. The authors 
also mention that the stripped gas tails behind the galaxy 
with a fraction of the galaxy's velocity and that the tails 
remain denser than the local ICM by a factor of a few. Our 
wakes behave differently. The stripped gas remains within 
~ 150 kpc from where it was stripped. Also the density in 
the tails of our galaxies are much lower, it exceeds 1.2 times 
the local ICM density only very close to the galaxy (see 
Fig.n^land lH]) . 

The different behaviour in the simulations of 
Ijachvm et al.l (|2007l ) has several reasons: The ram pressure 
interaction is described differently, their effective resolution 
in the tails may be lower than in our simulations, most 
of their clusters are more compact and have a lower ICM 
density, and they restricted the ICM-ISM interaction to the 
inner 300 kpc of the cluster. 

4.4 Comparison with observations 

4.4-1 Long HI tails in cluster outskirts 

IChung et all (|2007l l have searched for HI tails in the Virgo 
cluster. Of the ~ 50 targeted spiral galaxies, 7 revealed long 
one-sided HI tails. Surprisingly, these galaxies are at pro- 
jected cluster-centric distances of 0.6 to 1 Mpc. For their 
sensitivity limit of ~ 2 ■ 10^^ cm~^, the tail lengths are 
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~ 30 to 40kpc, while tail widths are similar to the cross- 
section of the remaining gas disc with respect to tail di- 
rect ion^jrh«_cohirnn_d« decreases gradually along the 
tail. IChung et al] l|2007l ) estimated likely ram pressures and 
gravitational restoring forces for the galaxies. Comparing 
the two forces, they found that 5 of these 7 galaxies could 
indeed suffer RPS. 

In our simulations, we have not modelled the Virgo clus- 
ter directly, but our cluster CI is also a compact cluster. In 
the literature, different parameters for the ICM distribution 
in the Virgo cluster are given (see Fig.[T}. For cluster-centric 
radii of 0.5 to 1 Mpc, our cluster CI is a factor o f abo ut 3.6 
denser than the version used in IVoUmer et all (l200lll. but 
agrees well with the parameters given by Matsumoto et all 
hood). Thus, we consi der a comparison of the tails observed 
bv lChung etahl l|2007l ) with our simulated tails in cluster CI 
reasonable. 

Adopting a column density limit comparable to the sen- 
sitivity limit of the observations, our simulations are able to 
produce tails with similar characteristics regarding length, 
width and structure also at large distances to the cluster 
centre (e.g. first row of Fig. |4]and first two rows of Fig. |6]). 
However, not all simulated galaxies at large cluster-centric 
distances show long tails. In the first row of Fig. (5] the 
tails are rather short, because here the ram pressure is very 
small due to the small orbital velocity of this galaxy. Thus, 
we conclude that it is possible that the tails observed by 
IChung et all (2007) arc produced by ram pressure stripping 
although the galaxies are still a large distance from the clus- 
ter centre. However, these galaxies need to have a substantial 
velocity component, i.e. about 1000 kms"^, in the plane of 
the sky. No t all g alaxies at large cluster-centric distances 
observed bylCh ung ct al. (2007) show long tails. The reason 
can simply be that not all galaxies have a sufficient velocity 
component in the plane of the sky. We expect that in more 
extended clusters like Coma, long ram pressure induced tails 
can be found at even larger distances to the cluster centre. 

Our velocity plots (Figs. [T^ to I14p are the ana- 
logue to the p ositio n- velocity diagrams (PVD) presented by 
IChung et all (|2007l ). However, given the sensitivity limit of 
observations, only the densest features of our velocity dia- 
grams will b e observable. In the observations presented by 
IChung et al.l (2007 ) . the gas tails appear as a "hook" that 
is attached to one edge of the galaxy and shows a slight 
slope towards the galaxy's systemic velocity. This hook has 
a small velocity width, < 200kms~^, in some cases even 
less. If this was the whole velocity width of the tail, this 
would be surprising, as one would expect a velocity width 
like in the galactic disc (~ 400kms~^). 

In our simulations, we can observe a similar "hook"- 
feature in our w^i-plots, if we consider only densest parts, in 
cases where the galaxy is not moving close to face-on. Also 
here, the velocity width in the hook is only < 200kms~^, 
which agrees with the observations. In the simulations, the 
"hook" is not as long as in observations. This feature is 
generated by the fact that in these cases the gas loss hap- 
pens mainly at one side of the galaxy, namely the part of 
the leading edge where the galaxy rotates al ong with the 
wind. Also most galaxies with tails observed bv lChung et al.l 
are not moving close to face-on. However, in our sim- 
ulations, this hook feature is not as clear as in the observa- 
tions. Moreover, nearly all of the tailed galaxies observed by 



IChung et all (HqO^) show this hook, whereas in our simula- 
tions the hook appears only occasionally and along prefer- 
able line-of-sights. In general, in the PVDs of our simula- 
tions, the tail is much more diffuse than the observed ones. 
A reason for this difference could be that real RPS does not 
proceed as turbulent as in our simulations, but that viscosity 
leads to a smoother flow. 

We note that in the observed galaxies, in almost all 
cases the velocity gradient along the tail/hook is not only 
towards the cluster mean, but also towards the galaxy's sys- 
temic velocity. Moreover, some of the observed galaxies have 
a very small radial velocity with respect to the Virgo cluster 
mean, so that at maximum a weak velocity gradient along 
the tail could be expected, if the gradient is due to the accel- 
eration of the stripped gas away from the galaxy. We spec- 
ulate that the gradient in the tail does not simply show the 
acceleration of stripped gas away from the galaxy, but that 
the deceleration of the gas disc's rotational component along 
the tail plays a crucial role. 

4.4.2 The case of NGC 4388 

The Virgo spiral galaxy NGC 4388 seems to be an excellent 
example of a ram pressure stripped ga laxy. It is known to 
have a ~ 35 kpc long tail of ionised gas (lYoshida et al ] |200l 
\2004) . Additionally, lOosterloo fc van GorkomI (120051 ) have 
found a ~ 120 kpc long HI tail associated with this galaxy. 

NGC 4388 has a high ra dial velocity of 1400 kms"^ (see 
IVollmer fc Huchtmeierll2003l and references therein) with re- 
spect to the cluster mean. Given its long tail in projection, it 
also must have significant velocity component in the plane of 
the sky. A reasonable assumption is that the velocity com- 
ponent in the sky is comparable to the radial one, which 
would lead to a total velocity of ~ 2000 kms~^ in the Virgo 
cluster rest frame, which is rather high. However, assuming 
a much smaller component in the plane of the sky would 
increase the true length of the HI tail a lot. Our current 
assumption already results in a true length of the HI tail of 
~ 1.5 X 120 kpc = 180 kpc. 

A reasonable value for the local ICM density at the 
position of NGC 4388 is a few 10"^* gcm"^. With a likely 
velocity of about 2000 kms~^, it experiences a ram pressure 
of about 10"^^ ergcm"'^. According to our simulations, this 
is enough to be stripped heavily, so regarding the degree of 
stripping, simulations and observations agree. 

Regarding the combination of orbital velocity and lo- 
cal ICM density, in our simulations, there is no case di- 
rectly comparable to NGC 4388. The closest our simulations 
get to the velocity-density combination of NGC 4388 is the 
peri-centre passage of runs C1-LG-FST-. . . or at ~ 800 Myr 
of run C3-SM-FST-MF. In both cases, the ram pressure is 
about 10~^^ ergcm"^, but the orbital velocity is only slightly 
below 1500 kms~^. In both cases, the mass loss per orbital 
length is about 2.7 • lO^M0kpc-\ Given that NGC 4388 
is probably a factor of ~ 1.4 faster than the two mentioned 
examples, also the mass loss per orbital length should be a 
factor of ~ 1.4 smaller, which would lead to a less dense tail 
than the ones in the mentioned runs. As an opposing effect, 
the line-of-sight and the tail of NGC 4388 make an angle 
of about 45°, so that projected densities appear a factor of 
1.4 higher than they would if the line-of-sight and the tail 
were perpendicular to each other. Also the projected gas 



© 2007 RAS, MNRAS 000,[TH23] 



18 E. Roediger and M. Briiggen 



mass per orbital length is about 1.4 times higher. As these 
two effects cancel each other, we expect column densities 
and gas masses per orbital length comparable to the ones 
in our simulations. Thus, our simulations predict that the 
column density in the tail near the galaxy should be near 
~ 10'^° cm~^, and decrease gradually with increasing dis- 
tance to the galaxy to below ~ 10^^ cm~^ at a distance of 
~ 50 kpc or latest 100 kpc behind the galaxy. At a distance 
of ~ 100 kpc, the gas mass per orbital length should corre- 
spond to the galaxy's mass loss per orbital length, which we 
estimated above to be about 3 ■ 10^ Mq kpc~^. 

However, the tail of NGC 4388 looks different. It is 
densest at a distance of ~ 100 kpc behind the galaxy, 
here the column density reaches ~ 10^" cm~^. The spatial 
width of this dense patch is about 5 kpc, thus this dense 
patch corresponds to a mass per orbital length of about 
4- 10® Mq kpc~^. Towards the galaxy, the column density as 
well as the mass per orbital length decrease. In HI, the tail 
even s hows a gap near the galaxy. In this gap. lYoshida et al] 
observed a plume of ionised gas. However, the total 
mass of this plume is only 10^ Mq, which corresponds to 
a mass per orbital length to 3 • 10 Mq kpc~^ and is much 
less than we expected. Thus, the distribution of stripped gas 
along the tail differs from what we expect. 

The width of the tail shows some flaring. In our simula- 
tions, we do not observe tail flaring in the sensitivity limit of 
these simulations. In the simulations, the tails show flaring 
only at much lower column densities. 

The ionised part of the tail NGC 4388 shows a filamen- 
tary structure . On th e basis of a deep spectroscopic study, 
lYoshida et all (|2004h derive the Ha velocity field, which 
turns out to be quite complicated. The authors identify sev- 
eral kinematic groups among the filaments. However, they 
also consider significant turbulent motions possible. They 
argue that the kinematic structure and metallicity of the ion- 
isation are not well-explained by a minor-merger scenario, 
but can be explained naturally by ram pressure stripping 
and an additional starburst driven superwind in the galac- 
tic nucleus. According to our simulations, turbulent motions 
in the tail are very likely. However, the mass of stripped gas 

near the galaxy is much less than we expected. 

The HI tafl of NGC 4388 ()Oosterloo fc van GorkomI 

l2005f ) remains hard to fit with the simulations also with re- 
spect to the velocity information. The velocity width of the 
tail is rather small over the whole length. Additionally, ac- 
cording to our simulations, the acceleration of the stripped 
gas towards the ICM velocity is finished approximately at 
a distance of 100 kpc behind the galaxy. For NGC 4388, 
however, the velocity of at the end of the tail differs from 
the galaxy's systemic velocity only by 500 kms~^. According 
to our simulations, a value comparable to the galaxy's ra- 
dial velocity with respect to the cluster mean, 1400 km s~^, 
would be expected. 

In summary, our simulations have difficulties to explain 
the details of NGC 4388's tail, where the most severe points 
are: 

• The distribution of stripped gas along the tail differs 
between observations and simulations. The simulations pre- 
dict much more gas near the galaxy. The mass per orbital 
length at a distance of ~ 100 kpc behind the galaxy is ap- 
proximately at the expected level, however, here the gas is 



confined in a smaller volume than expected and thus shows 
higher column densities than expected. 

• According to our simulations, the stripped gas should 
be decelerated to the cluster mean velocity at a distance of 
~ 100 kpc. According to the HI data, for NGC 4388, how- 
ever, the stripped gas seems to have lost only 1/3 of the ex- 
pected velocity change. In contrast to this, the Ha tail shows 
a velocity gradient as expected from our simulations. Radial 
velocities at the end of the Ha tail have already reached 300 
to 500kms~^. Also the turbulent structure of the Ha veloc- 
ity field is in agreement with our simulations. We note that 
here the HI and the Ha data seem to be in disagreement. 

There are several reasons that could cause the mis- 
matches between the simulations and the observations. E.g. 
the ICM may not be smooth but inhomogeneous, and there 
may be ambient motions in the ICM, thus leading to a 
special ram pressure history of this galaxy. The Virgo spi- 
ral NGC 4522 is another example that suggests ambient 
motions in the ICM, as its degree of stripping implies a 
much stronger ram pressure than expected at its position 
in the cluster (.VoUm er et al. 2006?). Moreover, NGC 4388 
has an active nucleu s and a starburst driven superwind (see 
lYoshida et al.ll2004l and references therein) , which will infiu- 
ence the dynamics of the ICM-ISM interaction. Addition- 
ally, our simulations cannot predict the temperature of the 
stripped gas and are thus unable to predict in which wave- 
lengths the stripped gas is observable. However, the close 
vicinity of the galaxy has been searched for ionised and neu- 
tral gas and the amounts found are significantly less than 
what we expect. 



4.4-3 Ha and X-ray tails 

Up to date, only a few galaxies with Ha and X-ray tails are 
known. There are two irregular galaxies in the cluster A 1367 
that h ave Ha tails of 50 k pc and 75 kpc length and 8 kpc 
width (jCavazzi et al.|[200lh . The tails are straight an d show 
some filamentary structure along the tail. lYaei et all (|2007h 
have observed the Coma galaxy D 100 and found an Ha tail 
that is exceptionally straight and narrow (60 kpc long, 2 kpc 
wide). The Ha tail of NGC 4388 has been discussed above 
in c onnection with it s HI tail. 

ISun et al] l|2006l ) have presented a long and narrow X- 
ray tail (70 kpc long, 7 kpc wide) for the galaxy ESQ 137- 
001 in the Coma-like cluster A 3627. Additionallv. lSun et al.l 
(j2003) have found a 40 kpc x 4 kpc Ha tail for the same 
galaxy that coincides spatially with the X-ray tail. The mass 
of the X-ray tail is estimated to be 10^ Mq, the Ha tail 
provides another few 10* Mq. Both, the X-ray and the Ha 
tail are nearly straight, although some very slight oscillation 
can be detected. In both wavebands, the highest emission 
is found near the galaxy, but there are some more brighter 
patches further downstream. If the brightness also traces the 
local gas mass, this behaviour would be expected if these 
tails are caused by ram pressure stripping. Also the total 
gas mass of the tails matches the ram pressure picture. 

A characteristic feature of the observed Ha and X-ray 
tails is their narrowness and straightness. Given that our 
simulations cannot estimate the temperature of the stripped 
gas, we are unable to predict the brightness distribution of 
our tails in Ha or in X-rays. However, the projected density 
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maps presented in this paper generally show much broader 
tails and make such narrow and straight tails difficult to 
understand. Although our simulations provide an important 
step towards a realistic description of RPS by modelling the 
flight through a cluster, still important physics seems to be 
missing. 



5 SUMMARY 

We have performed 3D hydrodynamical simulations of ram 
pressure stripping of disc galaxies on orbits through galaxy 
clusters. Here we have presented a detailed study of the 
wakes associated with these galaxies. 

For a sensitivity limit typical for current HI observations 
(10^^ cm~^), we find that a typical tail length in projected 
gas density is 40kpc. In this sensitivity limit, the width of 
the tails is similar to the galaxy's width with respect to the 
tail direction. Our simulations show that the density and 
length of the tail is not determined by the current ram pres- 
sure and thus temporal mass loss rate alone, but also by 
the galaxy's velocity. These two dependencies can be sum- 
marised in the mass loss per orbital length, which is the 
crucial parameter that determines the tail length and den- 
sity. 

If we consider projected gas densities below the sensi- 
tivity limit, the tails show a flaring width, where the flaring 
ratio depends on the galaxy's cross-section with respect to 
its direction of motion. Large cross-sections lead to strong 
flaring, while for small cross-sections the flaring is weak. The 
flaring is caused by turbulence in the galactic wakes. 

The velocity in our wakes shows a significant turbulent 
component of a few lOOkms"^. The stripped gas is fully 
decelerated to the ICM rest frame at a distance of ~ 100 kpc 
behind the galaxy. 

We find that the stripped gas does not follow the galaxy 
for a long way but remains within ~ 150 kpc from the po- 
sition where it is lost. We have investigated in how far the 
deceleration of the galaxy by the drag force can be a heat- 
ing source for the ICM and found that no significant effect 
is expected. 

As our simulations neglect thermal conduction, cooling 
and viscosity, they are unable to predict the wavelength at 
which the galactic wakes should be observable. Nonetheless, 
we compared our simulations with observations of galactic 
wakes. 

We can reproduce tails of 40 kpc length at cluster- 
centric distances between 0.5 and 1 Mpc in a compact clus- 
ter similar to the Virg o cluster. Thus we c onclude that the 
galactic tails found bv lChung etHI (|2007l ) could indeed be 
formed by ram pressure stripping. However, our simulations 
have difficulties in reproducing the velocity information of 
these galaxies. The observations suggest that RPS proceeds 
less turbulent than in our simulations. Also, the structure 
of the wake of the Virgo spiral NGC 4388 is hard to ex- 
plain by our simulations, not only due to shortcomings of 
our simulations, but also due to intrinsic discrepancies in 
this galaxy's data. This HI "tail" may not be caused by RPS 
at all, or this galaxy may have experienced a complex ram 
pressure history. Given that our simulations cannot predict 
the correct gas temperatures, comparisons to the Ha and 
X-ray tails are especially difficult. However, our simulations 
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Figure A3. Same as Fig.Hbut for run Cl-LG-FST-MF. The Ihs 
column is for the resolution described in Sect. 12. Il the rhs column 
for a run where the resolution has been improved by a factor of 
2 everywhere. 



do not seem to be able to match the straightness and nar- 
rowness of these tails. Also these observations suggest that 
some physics that were neglected in our model, e.g. viscosity, 
cooling, and tidal interactions, may play an important role 
for the detailed structure of ram pressure stripping tails. 
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APPENDIX A: RESOLUTION 

In this section, we show slices in the orbital plane that dis- 
play the local gas density and the fiow structure (Fig. IA2p . 
projected gas density maps (Fig. IA3|I and velocity plots 
(Fig.im. These Figs, are for the run Cl-LG-FST-MF, but 
once for the spatial resolution as described in Sect. 12. Il and 
once for a run were the spatial resolution as been improved 
by a factor of 2 everywhere. Figure ESl compares the distri- 
bution of stripped gas along the orbit. 

In all plots, the basic structure is same for both resolu- 
tions. We conclude that our conclusions are not influenced 
by the resolution. 
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Figure Al. Same as Fig. 1151 but for run Cl-LG-FST-MF. Tlie two top rows arc for a resolution as described in Sect. ^2A\ in tlie two 
bottom rows the resolution is improved by a factor 2 everywhere. 
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s in Mpc 

Figure A5. Distribution of stripped gas along orbit (compare to 
Figs. [16] and [nil: for run Cl-LG-FST-MF, averaged over 300 kpc, 
for two different resolutions. Colours code the resolution (LR is 
low resolution, HR is resolution improved by a factor of 2 every- 
where). The thick lines are distribution as seen in the simulation. 
The thin coloured lines are the predictions based on the numerical 
result of the bound gas mass as a function of covered distance - 
however — shifted by 150 kpc along the orbit. The thin black line 
is the prediction based on the analytical estimate of the strip- 
ping mass (see paper I). The zero-point of the x-axis is shifted to 
peri-centre passage. 
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